Post-Burst Quasi-Periodic Oscillations from GRO J1744— 28 

and from the Rapid Burster 



Jefferson M. Kommers, Derek W. Fox, & Walter H. G. Lewin 

Department of Physics and Center for Space Research, 
Massachusetts Institute of Technology, Cambridge, MA 02139 

Robert E. Rut ledge 

Max-Planck-Institut fur Extraterrestrische Physik, D- 85740 G arching, Germany 

Jan van Paradijs 

University of Alabama in Huntsville, Huntsville, AL 35812 
and University of Amsterdam, Amsterdam, The Netherlands 

Chryssa Kouveliotou 

Universities Space Research Association, Huntsville, AL 35800 



ABSTRACT 

The repetitive X-ray bursts from the accretion-powered pulsar GRO J1744— 28 show 
similarities to the type II X-ray bursts from the Rapid Burster. Several authors (no- 
tably Lewin et al.) have suggested that the bursts from GRO J1744— 28 are type II 
bursts (which arise from the sudden release of gravitational potential energy). In this 
paper, we present another similarity between these sources. Rossi X-ray Timing Ex- 
plorer observations of GRO J1744— 28 show that at least 10 out of 94 bursts are followed 
by quasi-periodic oscillations (QPO) with frequencies of ~ 0.4 Hz. The period of the 
oscillations decreases over their ~ 30-80 s lifetime, and they occur during a spectrally 
hard "shoulder" (or "plateau") which follows the burst. In one case the QPO show a 
modulation envelope which resembles simple beating between two narrow-band oscilla- 
tions at ~ 0.325 and ~ 0.375 Hz. Using EXOSAT observations, Lubin et al. found 
QPO with frequencies of 0.039 to 0.056 Hz following 10 out of 95 type II bursts from 
the Rapid Burster. As in GRO J1744— 28 the period of these oscillations decreased over 
their ~ 100 s lifetime, and they occurred only during spectrally hard "humps" in the 
persistent emission. Even though the QPO frequencies differ by a factor of ~ 10, we 
believe that this is further evidence that a similar accretion disk instability is responsible 
for the type II bursts from these two sources. 
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1. Introduction 

GRO J1744— 28 is a recently discovered accretion- 
powered pulsar which shows repetitive X-ray bursts 



( Kcpuveliotou et al. 1996~Strickman et al. 1996~Gilcs 



before the bursts. Subsequent RXTE observations 
showed that the bursts were sometimes preceded by 
steadily increasing variability in the persistent emis- 
sion, including "micro" and "mini" bur sts (piles & 



1. 1996; and references therein). Its unusual burst- 



ing behavior has been compared (see below) to that of 
the Rapid Burster (MXB 1730-335), which was dis- 
covered 20 years ago by Lewin et al. (1976). Analysis 
of the X-ray bursts from the Rapid Burster showed 
that it produced two distinct types of burst: type I, 
attributed to thermonuclear flashes on the surface of a 
neutron star; and type II, attributed to the release of 
gravitational potential energy due to spasmodic ac- 
cretion. The mechanism responsible for the type II 
bursts has not been fully understood, although it is 
almost certainly related to an accretion disk insta- 
bility (for reviews see Lewin, van Paradijs, & Taam 



1996Q . 



Strohmayer 1996; Lewin et al. 1996; Giles et al 



Extensive reviews of the Rapid Burster and its 
complex behavior can be found in Lewin, van Paradijs, 
& Taam (1993, 1995); see also references therein. 
Here we discuss the features of the type II bursts from 
this source that are relevant to our comparison with 
GRO J 1744- 28. The time intervals between type II 
bursts from the Rapid Burster range from ~ 10 s to 
~ 1 hr, with the shorter intervals being more com- 
mon. Burst durations range from ~ 2 s to ~ 680 
s. The burst repetition pattern is that of a relaxation 
oscillator: the fluence of a type II burst is roughly pro 



199.3,1995). Rapidly lepelilive type II buisLs had pie- 



viously been ubsei ved only fium the Rapid Bmslei , su 

if the bursts from GRO J1744-28 are also of type II — 
as convincingly argued by Lewin et al. (1996) — then 
a comparison of these sources may constrain theories 
of fee burst mechanism. 



porti onal to the time to the next burst ( Lewin et al 
1976; Lewin et al. 1995). Persistent X-ray emission 



between the type II bursts is observed following long 
(duration > 30 s) bursts. The persistent flux emerges 
gradually after high-fluence bursts and decreases prior 
to the next burst (Marshall et al. 1979 ; Van Paradijs 



The pulse period of GRO J1744-28 is 467 ms and 
the neutron star is in a 11.8 day binary orbit about 
a low-mass donor star (Finger et al. 1996). During 
the first 12 hours of bursting observed with BATSE 
on December 2, 1995, bursts occurred every 3 to 8 
minutes, and for one three hour period the burst in- 
tervals cl ustered around 172 ± 15 s ( Kouveliotou et 
al. 1996| ). Subsequently the burst intervals became 



Cominsky, fc Lewin 1979; Btella et al. 1988a). These 



longer and more erratic and the burst rate settled at 
about 30-40 p er day (corrected for Earth occultation 
and live time; 
were initially < 



Fishman et al. 1996 ). Burst durations 
-> 20-30 s, and settled down to ~ 5-10 s. 
Unlike the type II bursts from the Rapid Burster, no 
relationship between the burst fluence and the time 
to the next (or previous) burst has been reported^ 
(Kouveliotou et al. 1996; Strickman et al. 1996). 



Using the first observations of GRO J1744— 28 with 
the Rossi X-ray Timing Explorer (RXTE), Swank 
(1996) noted that the bursts were followed by a char- 
acteristic "dip" in the persistent emission level which 
took a few minutes to recover. No such dips were seen 



A preliminary analysis of the BATSE catalog of over 3000 
bursts from GRO J1744— 28 shows a statistically significant 
but weak correlation between the burst fluence scaled to the 
persistent emission level and the time to the next burst (Kom- 
mers et al. 1997). 



pre- and post-burst features are referred to as "dips" 
in the persistent emission. The spectrum of the per- 
sistent emission is relatively soft during the dip just 
after a burst. It then rapidly increases in hardness, 
remaining hard for ~ 1-2 minutes before gradually 
decreasing to again become very soft during the dip 
preceding the next burst ( )5tclla et al. 19881 ). The 
1-2 minute period of spectrally hard emission corre- 
sponds to a "hump" in the persistent emission light 
curve. 

Lubin et al. (1992) found "naked eye" quasi- 
periodic oscillations (QPO) following 10 of 95 type 
II bursts from the Rapid Burster observed with EX- 
OSATva August 1985. The oscillations occurred only 
during the spectrally hard humps (which immediately 
followed the post-burst dips). The frequency of the 
oscillations ranged from 0.039 to 0.056 Hz, with a pe- 
riod decrease of 30-50% observed over the ~ 100 s 
lifetime of the oscillations. The fractional root-mean- 
square (rms) variation in the oscillations was 5 to 
15%. In 8 of the 10 cases, the "naked eye" oscilla- 
tions were accompanied by ~ 4 Hz QPO with rms 
variations of 6 to 19% ( Lubin et al. 1992 ). 

Several authors have already noted similarities be- 
tween the bursts from GRO J1744-28 and the Rapid 
Burster. After the discovery of the bursts from 
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GRO J1744-28, Kouveliotou et al. (1996) suggested 
that the release of gravitational potential energy fol- 
lowing an accretion instability might be responsible 
for the bursts. Lewin et al. (1996) made a detailed 
comparison between the two systems. Noting that 
both are transient low-mass X-ray binaries in which 
the accretor is a neutron star, they concluded that the 
bursts from GRO J 1744- 28 must be type II based 
on (i) the hardness of the burst spectra, (ii) the lack 
of spectral evolution during the bursts, (iii) the fact 
that the ratio of integrated energy in the persistent 
emission to that in the bursts was initially too small 
to allow for a thermonuclear burst mechanism, and 
(iv) the presence of dips in the per sistent emission 
following bursts flLcwin et al. 19961 ). Subsequently, 
Sturner & Dermer (1996) reached the same conclu- 
sion. The possibility of thermonuclear burning in 
GRO J1744-28 has been discussed by Bildsten & 
Brown (1997). 

In this paper, we present another similarity be- 
tween the Rapid Burster and GRO J 1744- 28. Both 
sources show transient QPO during spectrally hard 
emission intervals following bursts. 

2. Observations 

Since January 18, 1996, RXTE has performed nu- 
merous observations of GRO J1744— 28. The data 
discussed here were taken by the PCA instrument and 
are publicly available. Between January 18 and April 
26, 1996, t here were 94 main bursts observed ( Giles 
2t al. 1996| ). To refer to these bursts individually, we 



number them sequentially from 1 to 94. 

The light curves of the bursts reveal a rich phe- 
nomenology in the burst profiles. Figure |l| (a) shows 
the "dip" in persistent emission which follows some 
bursts, as first noted by Swank (1996). A broad 
"shoulder" (or "plateau" ) of emission above the mean 
pre-burst level immediately follows some bursts, as 
shown in Figure @ (b). The shoulder occurs before 
the dip whenever both features are present. Some 
bursts show large-amplitude oscillations with frequen- 
cies ~ 0.4 Hz during their shoulder, as shown in Fig- 
ure [j] (c) ( Kommcrs et al. 1996| ). 

Oscillations can be seen during the shoulders of the 
light curves immediately following at least 10 of the 94 
bursts. These ten bursts were numbers 8, 12, 14, 18, 
43, 53, 65, 68, 77, and 94. Typically 5 to 15 cycles of 
the QPO are apparent, with as many as ~ 25 cycles 
seen in the case of burst 65. The low frequency of 



the oscillations and the low numbers of cycles make 
it difficult to use Fourier power spectra to study these 
QPO. We have instead found the mean periods of the 
oscillations by estimating the time intervals between 
successive maxima in the count rates. 

The mean frequency of the oscillations over the 
ensemble of bursts is 0.38 ± 0.04 Hz. The mean 
frequency of the oscillations after individual bursts 
varies from 0.35 Hz ± 0.02 Hz in burst 65 to 0.49 
± 0.03 Hz in burst 12. (The uncertainty in these 
figures represents the uncertainty in the mean fre- 
quency.) The period of the QPO following a given 
burst typically wanders non-monotonically about a 
mean period by ~ 20 percent over the lifetime of the 
oscillations. Using the \ 2 statistic we exclude (at the 
2a level or better in each of the 10 bursts) the null 
hypothesis that the period between successive QPO 
pulses is constant. 

To identify any overall tendency towards increas- 
ing or decreasing QPO periods, we looked for bursts 
where the Spearman rank-order correlation coefficient 
(r s ) indicated that the time intervals between QPO 
pulses were correlated (or anti-correlated) with arrival 
time. In no case did we find a positive r s , which would 
have indicated an overall trend towards increasing pe- 
riods. For bursts 18, 53, 65 (see Figure 2), and 94, we 
found an anti-correlation at the 2a level or greater. 
In these bursts the QPO period tends to decrease by 
10 to 20 percent over the lifetimes (~ 50 to 80 s) of 
the oscillations. 

The fractional rms variations of the oscillations 
range from ~ 5 to 13%. Using the counts to flux 
conversion 4.1 x 10~ 12 erg cm~ 2 s _1 per count (to- 
tal PCA band; |Gilcs et al. 1996 ) the rms variations 



in flux units range from 1.2 x 10 erg cm s to 
1.4 x 10~ 8 erg cm" 2 s" 1 . 

The shoulders during which the post-burst oscilla- 
tions occur are spectrally harder than the persistent 
emission. Figure || (top panel) shows a hardness ratio 
for burst 65. The hardness is defined as the count rate 
(averaged over 1.0 s) in the 3-5.6 keV range divided by 
that in the 2-3 keV range. The middle panel shows 
the power spectrum as a function of time. Higher 
power levels are shown with darker shades. The pres- 
ence of the QPO is shown by a dark band at 0.4 Hz 
lasting from the end of the burst until ~ 100 s after 
the peak. The bottom panel shows the burst intensity 
profile. The spectrum of the burst counts is clearly 
harder than that of the persistent emission, but even 
after the main part of the burst subsides the hardness 
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ratios during the presence of the QPO (20 to 100 s af- 
ter the burst) exceed those of the pre-burst persistent 
emission. 

The QPO following burst 65 show a modulation en- 
velope reminiscent of "beating" between oscillations 
at separate frequencies; see the bottom panel of Fig- 
ure ^. The count rates shown in the inset have been 
smoothed with a boxcar average to highlight the enve- 
lope of the (roughly) 0.05 Hz "beats". If this modula- 
tion envelope can be attributed to simple beating, the 
two narrow-band oscillations must occur concurrently 
and have roughly comparable amplitudes. An esti- 
mate of the two frequencies involved can be obtained 
from the mean frequency f ave = (/j + /2V2 and the 
beat frequency fbeat = \ fx ~ fs\- The mean frequency 
for the ~ 25 cycles of QPO following this burst is 
0.35 ± 0.02 Hz. The two frequencies which are beat- 
ing must then be roughly f x = f ave - fbeat/ 2 = 0.325 
Hz and / 2 = fave + fbeat/2 = 0.375 Hz. The aver- 
age frequency wanders by about 20%, so if beating 
is responsible for the apparent modulation envelope 
then the frequencies of the two narrow-band oscilla- 
tions wander as well. The presence of a similar phe- 
nomenon following other bursts is difficult to deter- 
mine because fewer cycles are available. 

Over the course of the January — May 1996 RXTE 
observations, the burst fluence, peak flux, and per- 
sistent emission level each decreased (approximately 
linearly) by a factor of ~ 4-5 (Giles et al. 1996). 



The average frequency and rms amplitude of the post- 
burst QPO show no significant (> 2a) correlations 
with any of these quantities. QPO with nearly the 
same frequency are seen in bursts for which the per- 
sistent emission level differs by a factor of ~ 4. 

Although our analysis focused on the December 
1995-May 1996 outburst of GRO J1744-28, we note 
that post-burst oscillations were also detected fol- 
lowing a burst in June 1996. This burst occurred 
when the source temporarily resumed bursting activ- 
ity at a lower level than the main December 1995- 



May 1996 outburst (Jahoda, Strohmayer, & Corbet 
19S(6|). Kommers et al. (1996) reported 0.4 Hz os- 
cillations during a ~ 25 s "shoulder" following the 
sixth of 7 bursts observed with the PCA on June 4, 
1996. The fractional rms amplitude of these oscilla- 
tions was 25 ± 5 percent. A second large outburst 
from GRO J1744— 28 has been in progress since De- 
cember 2, 1996, but we have not yet analyzed these 



observations ( Kouvcliotou et al. 
19S0). 



1997; Stark et al 



3. Discussion 

The ~ 0.4 Hz oscillations following some bursts 
from GRO J1744— 28 are reminiscent of oscillations 
that have been seen following some type II bursts from 
the Rapid Burster. This likeness between the two 
sources complements the comparisons made by Lewin 
et al. (1996). Although it is not certain that the post- 
burst oscillations are the same phenomenon in both 
sources, the following similarities are worth consid- 
eration, (i) The post-burst oscillations follow bursts 
that do not show the profiles or spectral evolution 
characteristic of Type I (thermonuclear flash) bursts, 
(ii) The oscillations occur during a period of spec- 
trally hard emission, (iii) When an overall change in 
QPO period is observed, it is a decrease: 10-20% for 
GRO J1744-28, and 30-50% for the Rapid Burster, 
(iv) The fractional rms amplitude of the oscillations 
is roughly 5-15%. 

There are some differences, however, (i) The fre- 
quency of the post-burst oscillations in GRO J1744— 28 
is roughly 0.4 Hz while in the Rapid Burster, the 
post-burst oscillations have frequencies of 0.04 Hz and 
are in some cases accompanied by 4 Hz QPO. (ii) In 
GRO J1744— 28, the oscillations occur before the dip 
in the persistent emission; but in the Rapid Burster, 
they appear after the dip flLubin et al. 1992 ). 

Cannizzo (1996a, 1996b) has shown that a sim- 
ple numerical model of a thermal-viscous accretion 
disk instability can reproduce some features of the 
GRO J1744-28 bursts. In his model, interplay be- 
tween radial and vertical energy transport in the disk 
causes oscillations in the ratio of gas pressure to total 
pressure (denoted by (3). These oscillations increase 
in amplitude and eventually become non-linear, lead- 
ing to a Lightman-Eardley instability (Lightman & 
Eardley 1974 ). For specific choices of the inner and 
outer disk radii (ri nner = 10 7 5 cm, r outer = 10 9 cm) 
and the form of the viscous stress, integration of the 
time-dependent model reproduces the ~ 1000 s re- 
currence times and the ~ 10 s durations of the bursts 
as well as the post-burst dips. Cannizzo (1996b) also 
points out that the increased variability before bursts 
may be related to the oscillations in (3 that precede 
the instability. These oscillations may occur in a vari- 
ety of modes, but in his idealized case appear to have 
frequencies of about 0.05 Hz (Cannizzo 1996b). It re- 
mains to be seen whether this model can explain the 
oscillations reported here. 

Abramowicz, Chen, & Taam (1995) have shown 
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that low-frequency oscillations in mass flow can arise 
in accretion disk-corona systems. In their model, 
the development of a Lightman-Eardley instability is 
moderated by energy dissipation in the coronal re- 
gion, which lies above the thin disk. They propose 
that the strong ~ 0.042 Hz oscillations seen in the 
Rapid Burster may arise from such a mechanism. The 
frequency of these QPO may be relatively insensitive 
to the accretion rate, since it is determined by several 
competing factors including the radius of the inner 
accretion disk. The spectrum of the QPO is expected 
to be hard whenever the QPO originate from the hot 
inner region (Abramowicz et al. 1995). 

The common feature responsible for oscillatory be- 
havior in these models is the presence of a mechanism 
which couples the vertical energy transport to the ra- 
dial energy flux. In Cannizzo's model, this mechanism 
also acts to produce the bursts. To make a quantita- 
tive comparison of these models with the characteris- 
tics of the post-burst oscillations in GRO J1744— 28 
and the Rapid Burster, one would have to know what 
QPO frequencies are expected based on the physi- 
cal parameters that distinguish the two sources. The 
magnetic field strength and inner disk radius in par- 
ticular should be quite different ( Lcwin ct al. 199(| ; 
Finger et al. 1996; ^turner fc Dcrmcr 1996| ), which 
might account for the factor of ~ 10 difference in 
post-burst QPO frequency. Detailed information of 
this kind has not yet been presented for the models 
of Cannizzo (1996a,b) or Abramowicz et al. (1995). 

The presence of transient post-burst QPO is an- 
other similarity between the bursts from GRO J1744— 28 
and the type II bursts from the Rapid Burster. Com- 
bined with the comparison of Lewin et al. (1996), 
this similarity is further evidence that the bursts from 
GRO J1744-28 are of type II. If the same disk insta- 
bility is responsible for the type II bursts from each 
of these sources, then the post-burst oscillations pro- 
vide another observational benchmark for models of 
the burst mechanism as well as the QPO. 
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(a) Burst 33 
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Fig. 1. — Light curves of 3 bursts from 
GRO J1744— 28. The horizontal and vertical scales 
are each different to highlight the feature of interest. 
Panel (a) shows the "dip" following burst 33. Panel 

(b) shows a "shoulder" after burst 65; it occurs be- 
fore a slight "dip". QPO with a mean frequency of 
<~ 0.35 ± 0.02 Hz are present in the shoulder but are 
not apparent at this scale and time binning. Panel 

(c) shows an example of strong QPO with mean fre- 
quency 0.37 ± 0.03 Hz occurring during the shoulder 
after burst 68. The termination of the shoulder and 
the "dip" occur after the plotted time interval. 




Time (s) since 30 Mar 1996 17:04:44 UT 



Fig. 2. — Hardness, power spectrum, and intensity 
profile of burst 65. The top panel shows the ratio of 
the 3-5.6 keV count rate to the 2-3 keV count rate. 
(No instrumental dead time corrections have been ap- 
plied, which may affect the hardness during the burst. 
The hardness during the shoulder is not significantly 
affected by dead time.) The middle panel shows a 
power spectrum with higher powers represented by 
darker shades. The lower panel shows the intensity 
profile. Notice that the 0.35 ±0.02 Hz QPO occur 
after the burst, and the spectrum during their pres- 
ence is relatively hard. The inset shows the portion 
of the intensity profile during which the QPO were 
observed. The data in the inset have been smoothed 
with a boxcar average to highlight the apparent beat- 
ing between oscillations at ~ 0.325 Hz and <~ 0.375 
Hz. 
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